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ABSTRACT 



A bubble transport equation has been developed that 
the bubble distribution, as a function of position, 
velocity, and size, must satisfy. An analytical model for 
bubble transport in the upper ocean is chosen and solutions 
are developed for this model. For a surface source only, 
these solutions compare favorably to experimental data in 
the near-surface region. The depth of this region of 
agreement depends upon the circulation field chosen, but 
is of the order of 3 meters. 
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I. INTRODUCTION 



Studies concerning the propagation of sound in the 
ocean and under water sound scattering from the sea surface 
have shown a discrepancy between existing theories and ex- 
perimental results. Because air has a markedly different 
density and compressibility than sea water and because of 
the resonant characteristics of bubbles, the suspended air 
bubbles in sea water have a profound effect upon underAv^ater 
sound and can explain part of the discrepancies that have 
been encountered (Ref. 1 and 2). There are theories that 
predict the effects of air bubbles on absorption, scatter- 
ing and attenuation of sound as Avell as their effects on 
sound velocity (Ref. 3) , but there is no theory at present 
which predicts the distribution of air bubbles in the 
ocean. Such a theory would enhance existing theories and 
result in better predictions on the behavior of sound prop- 
agation in the upper ocean. 

A fe\'j experimental studies have been conducted on the 
distribution of bubbles in the ocean. Blanchard and Wood- 
cock were pioneers in this field, using a bottle to scoop 
up Avater from breaking Avaves and counting the number of 
bubbles of different sizes (Ref. 4). Later studies con- 
ducted at the U.S. Naval Postgraduate School, Monterey, used 
acoustic techniques to determine the bubble population, 
taking advantage of the fact resonant acoustical scattering 
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and absorption cross section of bubbles are o£ the order 

of a thousand times the geometrical cross section. These 

experiments indicated that bubbles of radius 60 microns or 

_ ^ 

greater have a distribution depth dependence of Z 

whereas smaller bubbles follow closer the exponental law 
- z /L 

e ' , where L is between 5 and 9 meters (Ref. 5). 

In this paper, development of a bubble transport 
equation has been outlined that the bubble distribution, as 
a function of position, velocity and size, must satisfy. 

An analytical model for bubble transport in the upper ocean 
is chosen and solutions are developed for this model. 
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II . TRANSPORT EQUATION 



A heuristic approach is used here to develop the 
bubble transport equation, since the results are the same 
as with a more rigorous approach and the mathematics is 
less involved (Ref. 12). Consider the distribution 
T (r , V , 1 , t)d"^rd vdl which defines the number of bubbles 
with radii in dl about 1 and with velocity in d"^v about v 
which are in a volume d r about r at time t. If collisions 
among the bubbles are neglected, then the follov/ing equa- 
tions can be written for the seven-dimensional space: 

1 = Ijiet flux in + S(f,v,l,t)J d^rd^vdl (1) 

7 7 

where S (r , v , 1 , t) d rd vdl is a distributed source of bubbles 
introduced in d^r about f at time t with radii in dl about 

7 

1 and velocity in d"^v about v. 

Referring to figure 1, which represents a volume ele- 
ment in real space, the net flux in, per unit volume, for 
the x-direction, is: 

fluxin-fluxout .. VY(x)TCx)dydz-V y(x+dx)T(x+dx)dydz 3 (v^T) 
volume dxdydz =‘9x 

Similiarly, the net flux in per unit volume for the y and 

j- 4.- • 3(''^u>i'') , 3(''^zT) T 

z-directions is - — r and - —5= 2. respectively. 

3y ■ 3z ^ 

Combining the above results, the total net flux in per unit 
volume for real space is - Vj.*(vT), where the subscript r 
on the del-operator means the operation is taking place in 
the real space. A similiar approach is used to obtain the 
total net flux in for velocity and size space: - (aT) 
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Figure 1 

Elemental Volume Sho^^^ing "Bubble" 
Streaming in Real Space (x,y,z). 
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for velocity space and - 9(R'l') for size space, where R=dl/dt. 

The total net flux in for the seven-dimensional space under 
consideration is the sum of the net flux in for real, 
velocity, and size spaces. Substituting this into equation 
(1) gives the general bubble transport equation: 



= 

9t 



_ V^. (aT) _ 9(RT) 

r V 



+ 



S 



( 2 ) 



The remainder of this paper will consider the special 
case of steady-state with no distributed sources, for which 
equation (2) can be written 



v.V^T + i.V Y + MI = - fv ,a + Y f3 

^ V 91 Lv aiJ '■ 

since f and v are independent variables. Although 
distributed sources will be neglected in this thesis, 
surface sources can be introduced as boundary conditions. 
It should be pointed out that the non-zero right hand 
side of equation (3) is peculiar to the bubble transport 
problem and arises becouse a is a function of v and R 
is a function of 1. 
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III. ANALYTICAL MODEL FOR TRANSPORT EQUATION 



Before solutions to equation (3) can be obtained, a 
realistic model for the bubble transport equation must be 
developed, with particular attention to formulating ex- 
pressions for the bubble acceleration (a), the time rate of 
change of the bubble radius (R) and the circulation field 
(V) . 



A, MODEL ASSUMPTIONS 

The model used in this development assumes Reynold's 
numbers less than one, which is the Stokes' Law Regime. 

With regard to this assumption two comments should be made. 
First, for small Reynold's numbers, bubbles will behave 
like solid spheres with little or no distortion and secondly, 
this places an upper limit of about 150 microns on the size 
of bubbles that can be considered. For Reynold's numbers 
much greater than one, (larger bubbles) the resistance of 
the turbulent wake v\?ill have to be included in the equation 
of motion, along with Stokes' drag force (Ref. 6). 

Surface tension has been neglected in the model, so the 
pressure inside the bubble is equal to the hydrostatic 
fluid pressure. Since the pressure from surface tension 
Zy/l, where y is approximately 74x10*2 NT/m for sea water, 
this assumption will introduce an error of less than 5 
percent for radii greater than 30 microns. Thus, where the 
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Reynold's number places an upper bound on the radius, neg- 
lecting surface tension places a lower bound on the bubble 
radius under consideration. 

The volume of a gas bubble in a liquid varies primarily 
through the effect of two influences. First, gas contin- 
uously leaks out of the bubble to dissolve in the surround- 
ing liquid. If the bubble was stationary, the ensuing loss 
of mass would necessarily lead to its ultimate disappearance. 
On the other hand, as a bubble assends or desends in the 
liquid, the ambient hydrostatic pressure decreases or in- 
creases and so does the pressure within the gas bubble, thus 
effecting the size. 

The effects of the interaction between gas diffusion 
and change in hydrostatic pressure can best be seen by 
looking at the ideal gas law: 




rTN/P 



where N is the total number of gas molecules, r is the gas 
constant, T is the temperature, and P is the pressure. Now, 
taking the time derivative of the ideal gas law equation 
and assuming isothermal changes. 



4tt1^R = ^ PdN/dt - Ndp/dt 
P 



and rearranging: 



R 



1 dN/dt 
3 N 



dP/dt 

P 
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The above development assumes a single gas inside the 
bubble vice some mixture of gases. In equation (4), the 
first term is a result of gas diffusion whereas the second 
term brings in the effects of compression. 

Due to the relative motion between the water and bubble, 
assume that the composition of the water in respect to dis- 
solved gas is everywhere uniform and saturated at atmospheric 
pressure, , except in a thin shell surrounding the bubble. 
Assume that the gas in the bubble is uniform right up to 
the gas-water interface, and that the liquid in contact 
with the bubble is saturated at the pressure, P, inside the 
bubble. Further assume that diffusion gradient is uniform 
throughout the thin shell. With these assumptions. Pick's 
First Law can be used, which relates the amount of material 
which diffuses through an unit area to the product of the 
concentration gradient and diffusion constant: 

dWdt = - 647t1^(P-P^) . 

In this expression 6 is a proportionality constant which 
depends on d, the thickness of the shell; on A, the dif- 
fusion constant of the gas in the water; and on a, the solu- 
bility of the gas in the water,, in accordance with the 
relation : 



6 = Aa/d. 



12 










# 



Using Pick's First Law with the hydrostatic pressure rela- 
tionship, equation (4) becomes: 

V 

R = 6rT — + - — (4) 

(a^-z) 3 (a^^-z) 

where P^/p^^g and is the bubble velocity in the 

z-direction. In reference 7, the expression 6rT was deter- 
mined experimentally to be approximately a constant, 

10”^ m/sec. 

B. EQUATION OF MOTION 

When air bubbles are entrained in a circulation field, 
there are three forces acting on them along with the bouyant 
force. Figure 2 illustrates an air bubble in a circulation 
field and shows graphicly the various forces which need 
to be considered in order to determine the bubble 
acceleration . 

If the bubbles were completely entrained by the fluid, 
i.e., bubble velocity = circulation field velocity, they 
would experience the same force as would the fluid enclosed 
in the same volume, dV/dt. However, complete entrain- 

ment can only occur if p = p^^, which is far from being true 
for air bubbles. Since the bubbles are only partially en- 
trained, the fluid flows past the bubbles creating a rela- 
tive velocity (u) and thus introduces a drag force. As 
mentioned earlier, small Reynold's numbers have been assumed 
and thus the drag acting on the bubble can be described by 
Stokes' Law, 
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Drag Force 

Adhering Mass 
Force 



Entrainment Force 



Figure 2 

Free Body Diagram of a Bubble, 

Showing the Forces Acting on 
the Bubble in a Circulation Field 
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^Drag 

The third force acting on the bubble is the innertial 
adhering mass force. When a spherical shape moves through 
the water, it drags along an amount of water equal to a 
fractional amount of the volume of the sphere. In the 
equations that follows, sigma (a) represents this fractional 
part, and in reference (8), sigma has been derived as 0.5 
for a sphere. 

Referring to figure 2 and using Newton's Second Law of 
Motion, the equation of motion for an air bubble in a cir- 
culation field can be written as: 

“ " Pw^di ■ "^^Pwir ■ ^ (P,,-P)^Sk* 

Using and dV/dt = (vV^)V yields 

a = - ^ u + £.k + 2 j!i£(v‘ V r") V (5) 

^2 o o 

where a = 912 nj Pp^^ and u=v-V. 

C. CIRCULATION FIELD 

Little is known about the circulation field in the 
upper ocean. Although some work is being conducted at pre- 
sent at the U.S. Naval Postgraduate School, Monterey. In 
reference (6), Levich develops a two-dimensional equation 
which describes the particle motion by wave action in an 
ideal liquid, and is the basic form used here to describe 
circulation field: 
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V = V^e 



COS kxk. 



+kz 



/s 

sin kxi 



V 

o 



C6) 



This equation gives a periodic form for the horizontal 
direction as well as an exponental decay in depth, which 
might be expected for the physical situation in the upper 
ocean. It is also well to note that the condition for an 
imcompressible fluid has been satisfied v\rith this equation, 
i . e . , V • V=0 . 

Figure 3 shows the cell structure implied by this cir- 
culation field. For this model, any two adjacent cells are 
identical to every other pair of adjacient cells that are 
separated from the first pair by an integer number of 
wavelengths. This implies that the boundary conditions must 
be periodic unless unphysical surface sources are to be 
assumed . 
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kx^ 




Figure 3 

Cell Structure o£ the Circulation 
Field Model 
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IV. SOLUTION TO TRANSPORT EQUATION 



In this section, a look is taken at the bubble transport 
equation with the model that was presented in the previous 
section. The main interest will focus at determining the 
bubble distribution as a function of depth and bubble size. 

Substituting into the bubble transport equation the 
expressions for a and R from the preceeding section results 
in the following equation: 

v-VrT + a-VyT + R9T/91 = [^‘^/l^ - ''z/3 (a^ -z)] T 



By defining 



(J)(Z,X,V^,V^,1) = 



T (r , V ,1) dydv 



y 



and using the two-dimensional circulation field, the 
above equation can be integrated over the y and Vy directions 
and the problem reduces to solving the partial differential 
equation : 



+ Vzii + axii .. ajii RS* = 



9''x 






2a 



V, 



^2 ■ 3Ta^-z) 



^ ( 7 ) 



zo ’ 



with the boundary condition <{)(0,x,v ,v ,1) = (0,x ,v 

X z o o 

Vxo»lo^ corresponding to a surface source. 

Equation (7) is a first order, linear partial differen- 
tial equation in 5 variables and can be solved by the method 
of characteristics (Ref. 9). Using this method, the partial 
differential equation is equivalent to the following 5 char- 
acteristic ordinary differential equations in terms of the 
parameter z, V'/hich was chosen for mathematical convience: 
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z = 0 


dx 


Vx 




3T " 


^z 


x=Xo 


II 


^x 

^z 


^x=^xo 




3z 




dz 


''z 




^ - 


R 


l = lo 


dz 


'^l 




d(J) 


2a -1 




dF = 


l^Vz 3(ai-z)_ 


(p (b=(t>n(^ jV ,V ,1 ) 

T ro ^ 0 ^ XO ZO 0 

' 



A close look at the characteristic ordinary differential 
equations reveals that they are simply the bubble dynamic 
equations . 

Solving these equations for <p v\;ould be straight forward, 
except the equation for v^ is non-linear and analytic solu- 
tions have been found only for very special cases (Ref. 10). 
Thus, approximate methods must be considered to obtain v^. 
The method used for this paper to solve for v^ was to 

assume v = V + u . Appendix A outlines this method with 
z z z 

the solution for u^ as: 



C 

= e” Jo 1 


J 0® lo 2 " 


g 


+ ifLc 






1 '^z 




^Wz 


9x 9z/ 



Since the drag acceleration on the bubble is large, the 
bubbles will be entrained by the fluid within the first few 
centimeters of travel and u^ thus approaches a terminal 
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velocity. Referring to the computer output for u^, shows 



that the assumption V >> u can be made and also u is 

independent of thus implying that is independent 

of V . 
zo 

By a similiar approach, it is shown in Appendix A that 
~ • Now, the characteristic equation involving x can 

be written as: 



This equation implies that the bubbles follo\v trajectories 
in the fluid and form streamlines as was indicated on figure 
3. Thus choosing a x^ one can describe the path which the 
bubbles follow. 

The solution of the characteristic equation for 1 is 
straight forward with 



The second term in this equation is a result of gas diffu- 
sion and limits the bubble lifetime and hence the depth 
that surface generated bubbles can penetrate. 

From the last characteristic equation, <}> takes the 
general form': 




and using the chosen circulation field. 



X 1/k sin ^ (e ’^^s inkxQ) . 
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where (}) as a function of C^>z,v v„,l) is obtained from the 
relationships between the characteristics and their variables 
Assume a surface source (J)^ that is separable, 

and for convience, let 



Since and are nearly independent of and this 

choice will not strongly affect the solution. Therefore, 



4.(x,z,v V ,1) = X (x )LC1 )6(v -V )6Cv -v ) 

X Z O O O XO X ZO z 



|(l-z/a^ 



1/3 

^EXP 



[^ 20 . 
Jol^v 



dz' 



Now, define the "bubble density" as a function of position 
and size. 



$(x,z,l) = (pCx,z,v ,v ,l)dv dv 
. i X z X z 

and integrate the expression for (}> over and v^: 

00 

$(x,z,l) = II ^Cx^)L(l^)6Cv^^-v^)5(v^^-v^) j^Cl-z/aiK' 



EXP 



f^ 2 a , , 

Jol^v 



dv^ dv, 

X Z J 



To take advantage of the 
integration are changed, 
the Jocobian: 



delta-functions, the variables of 
dv^dv^ = JdVj^Q ^"'^zo’ J is 



J 



^(^xo’^zo^ 
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The equation of interest now has the form: 

$(x,z,l) = X(x )L(1 ) (1-z/a.)^'^^ EXP f^ia ^z'J (8) 
0 0 ^ ■'OlZv^ 

Because of the number of mathematical steps involved in 

deriving the Jocobian, the development is found in Appendix 

B. The final approximate result is: 



V 



zo 



fZ 



V, 



2a 



where , 



°1^V2 



6 rt 



dz' EXP 






? = +2“ (v -Vj„^ 
l-JVz z z ax 



,^z' (l-z’/a^) 



dz 



o Vz 

Thus, the equation for the distribution of bubbles with one 
initial surface velocity can be written as (Ref. 12): 



i(x,z,v « X(Xo)LClo)(l-^/“l)^^^^ (9) 

Looking at equation (9), the term ef^^dz’ will be the pre- 
dominate feature in solving for $. Two things should be 
noted about First, 5 is seen to be proportional to 1/V^^ 

when the approximation is used. This result illustrates 

the importance of selecting a good circulation field model. 
Secondly, ^a6RT and shows that gas diffusions is very 
important in the solution for the bubble distribution. 



22 . 



V. RESULTS 



A computer program was written for calculating equation 
(9), to determine if the results were physically realistic. 
In this program, the assumptions were made that X(x^)=l and 
L (1 q) =K 1^ ’ (Ref. 11). For the latter assumption, K=1 
was used in the computions since the interest was in the 
shape of the distribution curve rather than the actual den- 
sities. Calculations for three circulation field veloci- 
ties were made and the following graphs have been plotted: 

^max • ^o 

$ (0,z,l) vs. Z for various radii 

$ (0,z,l) vs. 1 for various depths. 

Zmax maximum penetration depth for a given size 

surface generated bubble. All calculations were made for 
Xq= 0, which implies a straight line trajectory downward. 

If some other x had been chosen, line integration would 
had to be performed on the characteristic curves in 5-space 
(x,z,v^,V3,l) . 

The numerical technique used in computing the integrals 
was 10-point Gaussian Quadrature, which approximates the 
integrand by a Legendre polynomial series. Although this 
technique was fairly accurate for most of the integrals, it 
could not be used in computing the relative velocity (u^) . 
For this calculation, it is noted that the integrand is 
essentially zero except for the last few percent of the 
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interval (0,z) where it increases dramatically. was 

therefore evaluated by using double precision 32-point 
Gaussian Quadrature on the latter portion of (0,z). 



In computing the relative velocity, the approximation 
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ze 



-kz 



, x=0, 



was used which is accurate provided |Z|< 10 meters. Also, 
since the calculations involved taking the difference be- 
tween two such integrals, very little error was introduced. 

It is appropriate now to make a few comments about the 
results that have been plotted. Figure 4 is a graph of 
maximun depth versus inital bubble size. This graph shows 
that bubbles have a finite life time as a result of gas dif- 
fusion and that smaller bubbles introduced at the surface 
will not penetrate as deep as larger bubbles. 

Figures 5 through 7 are plots of (0,z,l) versus depth 
for various bubble radii. These graphs show the influence 
that the circulation field velocity has on how deep the 
bubbles will penetrate. It should also be noted that even 
though the surface density falls off as 1^'^*^, there is a 
certain depth where the larger bubbles become more populous 
than the smaller bubbles because of gas diffusion. 

Figures 8 through 11 are plots of $ (0,z,l) versus 
radius for various depths. The main feature about these 
graphs is that the initial distribution appears to persist 
to some depth which is a function of the circulation field 
velocity. 
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Figure 4 

Curves for Maximun Penetration vs.. 

Bubble Radius for Various Circulation 
Field Velocit ies 
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Figure 5 

Bubble Density vs Depth for Various 
Size Bubbles, v^=-0.5 m/sec. 
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Figure 6 

Bubble Density vs Depth for Various 
Size Bubbles, v^=-1.0 m/sec. 
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dl (number /m'^) with dl = l micron 




Figure 7 

Bubble Density vs. Depth for Various 
Size Bubbles, Vq=-1.5 m/sec. 
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Figure 8 



Bubble Density vs Bubble Radius 
for Various Depths, Vq=- 0.5 m/sec. 
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Figure 9 

Bubble Density vs Bubble Radius 
for Various Depths, v^=-1.0 m/sec. 
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Figure 10 

Bubble Density vs Bubble Radius 
for Various Depths, Vq=-1.5 m/sec. 
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VI . CONCLUSIONS 



There is little experimental data available to compare 
the results to, however, in reference (5) where measurements 
were taken acoustically, general shapes for bubble distri- 
butions are given. One should note that the experimental 
data in reference (5) is not for isoradius bubbles since 
the resonant frequency is a function of depth, however, a 
good comparision can be made within a few meters of the 
surface . 

Although the bubble distribution as computed does not 
follow a simple mathematical function, the results for the 
first few meters can be approximated by exponental laws, 
depending on the circulation velocity and bubble size. 

Looking first at figure 5, the distribution drops off 
rapidly with depth pointing to the fact that for small cir- 
culation velocities, a distributive source is needed since 
the depth penetration is small. For V^=-1.0 m/sec, the 
curves can be approximated for the first 3 meters by the ex- 
ponental laws, e’^^^ for small bubbles and e’^^^ for large 
bubbles. The curves for V^=-1.5 m/sec follow fairly close 
e~z/9 for the smaller bubbles when z>-3.5m and for 

larger bubbles when z>-4.5m. The curves for V^=-1.0 and 
VQ=-1.5m /sec compare favorably with the experimental data 
mentioned, especially for the smaller bubbles. 

Even though there is fairly good agreement in the 
distribution near the surface, a definite discrepancy does 
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occur at the deeper depths. Probably, the major reason for 
this discrepancy is the fact that no distributive or bottom 
sources were considered in the calculations. The data in 
reference (5) was taken in shallow water and calm seas, so 
inclusion of such sources into the model should improve the 
agreement between the solution and the experimental data. 
Another possible reason for the discrepancy is that a large 
portion of the Jacobian was assumed small and neglected. 

This term was not looked at in detail and could be of some 
importance at the greater depths. 

Since the results do agree in part with experimental 
data, the approach taken in this paper does appear feasible 
for determining the bubble distribution in the upper ocean. 
For further studies in this area, it is suggested that solu- 
tions to the bubble transport equation be extended to include 
calculations for x?^0, models with distributed sources and 
surface tension and better models for the circulation field 
and gas diffusion. 
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APPENDIX A 



Solutions for the Relative Velocities 



In this appendix, solutions are obtained for the rela- 
tive velocities from the characteristic equations. Consider 
first the characteristic equation for v^ , which is non-linear 
Write differentiate: 

^^7 = dVz ^ duz 
dz dz dz 

Substituting this result into the characteristic equation 
for v^ gives 



du. 



g ^z _ g 
1 2 V 7 ~ ovr 



V, 



V. 



A formal solution to this equation can be 

•z ,7 . .dz". 



fz— dz ' 
“z = ^ 



P ' r^v. 

Jo 



r 






1 

+ a 



L 



3Vz 3Vz 
written : 







3Vz\ 




3x 


+ 3z / 
J 



dz' + uzo > 
or, rewritting. 




For computational purposes, the’ approximation v = V will 
be made where appropriate, since the above expression is an 
implicit equation. This approximation turns out to be valid 
as indicated in the Results Section. 
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In order to solve the characteristic equation for v , 
the same approach is used as described above, with the 
result : 



A close look at this equation reveals that u^ is approximately 
zero since the exponental term shows a sharp decaying 
feature and the circulation field chosen does not c’ ige 
rapidly. Thus the approximation can be used i Che 

calculations . 
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APPENDIX B 



Development of the Jacobian 



This section is devoted to the mathematical manipula- 
tions involved in developing an expression for the Jacobian 
required for equation ( 8 ) . The follov^^ing are the functional 
relationships among the variables and characteristics that 
need to be considered: 



V. 



“ v^Cz.Xo.Vj^o.v^o.Iq) 



V = V (z ,x ,v ,v ,1 ) 
z z ^ ’ o xo zo o 

1 1 ("’^O’^XO’^ZO’^0^ 

X = X (z,x^,v^q,v^q,1q) 

^X ~ ^x ’^O ’^XO >^ZO ’ ^o »^xo^ 

v_ = V (z,x V ,v ,1 ,V ) . 
z z o xo zo o zo 

For mathematical convience, define: 

= 8 Vj^/ 8 vxo ; Jz - 9Vx/3vzo ; J 3 = 

J 4 = dV2,/^^zO ’ *^5 - ^^/^^XO ’ '^6 ~ ^^/^^ZO* 

Thus , 

, _ aCvx.vz) 

° 3(''xo>Vzo) ' 1 J 4 - J2J3- 

Now, using the characteristic equations for v ,v and x 

X z 

and noting that, for example: 



the following equations can be written: 
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For ease of manipulation, define: 
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= 2a Cvx-Vx) ; q2 ^ 2a (v^-Vz) 

l^vz - 

The six equations can now be written as: 



dJ 



dz^ AJ]^ + BJ3 + E Jr cj, — L 

^ J- 3 v 
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Taking the moments of the first four equations 
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Jl "4 . CJ 1 J 4 ^ DJjJj ^ FJjJj ^ qj3i Jj. 

z o 



Now, combining the above equations give the result: 
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where 

Q E eC'^4'^5 - - '^2*^5^ ' 

Using the equation for 1 that was developed earlier and taking 
the partial derivatives gives: 
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If gas diffusion had been neglected in the development, the 
above two equations would have been zero and the solution to 
the transport equation would become unphysical. 

Substituting the above results into the differential 
equation for J results in: 
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where now, 
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For simplification and as a first approximation assume 
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and the solution to this differential equation is: 
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RELZ AS A FUNCTION OF Z(METERS) AND L(MICRONS) 
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